Spin-orbit coupling induced anisotropy effects in bimetallic antiferromagnets: A 
route towards antiferromagnetic spintronics 
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Magnetic anisotropy phenomena in bimetallic antiferromagnets Mn2Au and Mnlr are studied 
by first-principles density functional theory calculations. We find strong and lattice-parameter 
dependent magnetic anisotropics of the ground state energy, chemical potential, and density of 
states, and attribute these anisotropies to combined effects of large moment on the Mn 3d shell 
and large spin-orbit coupling on the 5d shell of the noble metal. Large magnitudes of the proposed 
effects can open a route towards spintronics in compensated antiferromagnets without involving 
ferromagnetic elements. 
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The introduction in 1991 of a hard-drive read-head 
based on the ferromagnetic anisotropic magnetoresis- 
tance (AMR) effect represented a major step in the 
field which later on became known as spintronics [![. 
The basic physics of the AMR, namely the bulk na- 
ture of this magneto-transport effect and its subtle 
spin-orbit coupling (SOC) origin, have limited the 
magnitude of the effect and the potential of scaling 
the AMR devices into small dimensions required for 
current high-density magnetic storage applications. 
These limitations were partly overcome by the dis- 
covery of the giant and tunneling magnetoresistance 
(GMR and TMR) effects which are interface trans- 
port phenomena and which rely primarily on the more 
robust magnetic exchange interaction. To date, stud- 
ies of magnetoresistive effects potentially suitable for 
AFM spintronics have focused on AFM counterparts 
to the GMR While the viability of this ap- 

proach is yet to be discerned it has been acknowl- 
edged that the requirements on the structural quality 
and the coherence of transport through interfaces in 
the GMR devices are significantly more stringent in 
the case of AFMs. 

This paper aims to open an alternative route to- 
wards AFM spintronics which reintroduces the lead- 
ing role of SOC. In this approach, the stringent re- 
quirements on the GMR/TMR are circumvented by 
considering instead the tunneling AMR (TAMR) and 
the Coulomb-blockade AMR (CBAMR) effects. We 
also demonstrate that SOC can be employed to con- 
trol magnetic anisotropies in the AFM in a way that 
leads to reorientation of the staggered moments re- 
quired for observing the above anisotropic transport 
effects. 

Previous studies of the TAMR and 



CBAMR [H E| in ferromagnets have shown that 
anisotropic magnetoresistance phenomena can be ex- 
tended from bulk to nanoscale devices, can have large 
magnitudes, and do not require spin-coherent trans- 
port throughout the structure. In transition metal 
ferromagnets a generic principle has been outlined, 
based on studies of magneto-crystalline anisotropies 
and of the TAMR and CBAMR @, M, G3], that 
the magnetic anisotropy phenomena are maximized 
in bimetallic systems combining large spontaneous 
moments on the 3d shell of a transition metal and 
large magnetic susceptibility and SOC on the 5<i-shell 
of a noble metal. Since Mn carries the largest moment 
among transition metals and most of the bimetallic 
alloys containing Mn order antifcrromagnetically, the 
goals of strong magnetic anisotropy phenomena and of 
AFM spintronics appear to merge naturally together. 
In our relativistic ah initio study we consider the 
Mn2Au AFM for which recent theoretical calculations 
predicted [13| record Neel temperature (> 1500 K) 
among Mn-based AFM alloys. The generic nature 
of the proposed anisotropy phenomena is confirmed 
by calculations in the conventional bimetallic AFM 
Mnlr. 

We begin with a detail analysis of magneto- 
crystalline anisotropies in M^Au. Calculations are 
performed considering the MoSi2-type bet structure 
with experimental lattice parameters, a — 6.291 a.u. 
and c = 16.142 a.u. The tetragonal unit cell with 
two formula units (f.u.) and an AFM arrangement of 
magnetic moments are shown in Fig. 1. This staggered 
magnetic structure was obtained from the theoretical 
study [HI of the magnetic interactions, and leads to 
the AFM stacking of the Mn layers with no frustra- 
tion. 
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FIG. 1: Total DOS (per f.u.), spin-resolved Au-atom and 
Mn-atom projected DOSs, and a schematic view of the 
crystal and AFM structure of Mn2Au. (Mn-I denotes one 
of two Mn-sublattices with opposite spin-polarizations.) 



We employ the relativistic version of the full- 
potential linearized augmented-planewave (FP- 
LAPW) method in the local spin-density approxima- 
tion. The SO coupling is included in a self-consistent 
second- variational procedure fl4j . This approach is 
typically very accurate for itinerant metallic systems. 
The calculated total and atom resolved DOSs for mo- 
ments aligned with the [001] axis are shown in Fig. 1. 
We find that Mn atoms carry spin Ms = ±3.2 \ib 
and orbital Ml = ^0.013 /is magnetic moments and 
there are no magnetic moments on the Au atoms. 
When staggered magnetization is aligned along the 
[100]/[010]-axis there is no change in Ms of Mn and 
the magnitude of \Ml\ = 0.007 \ib is reduced. 
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TABLE I: Relativistic full-potential DFT calculations of 
the uniaxial Ki, and forth-order out-of-plane, and in-plane 



Km cos 40) sin 2 9]. Here K 2 ± is the uniaxial MAE 
constant, and K 4 ± and Km\ are the forth-order out-of- 
plane and in-plane MAE constants, respectively. The 
advantage of this approach is that it allows us to split 
the total MAE into the element-specific contributions 
from different atoms in the unit cell. The total and 
clement specific anisotropy constants are shown in 
Tab. I. They are calculated with accuracy better than 
0.01 meV/f.u. 

The largest uniaxial anisotropy constant, K 2 a_, is 
strong and clearly dominated by contributions from 
the Au sublattice. The leading role in the MAE of 
Au-atoms which carry zero net moment may appear 
counterintuitive and requires more detail inspection. 
In the AFM state the Au orbitals point towards the 
neighboring Mn moments with opposite orientations 
and strong hybridization between Mn and Au creates 
oppositely polarized parts of the Au orbitals. Con- 
sidering separately exchange fields produced by the 
two Mn sublattices, a =t>b with opposite strongly 
localized moments ( ±M 3d ), each of the sublattices 
induces magnetic moment on the itinerant 5d 



electrons of Au atoms. We can write 12, 15 
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where we sum over the Mn atoms in the sublattice a, 
^3d-5d 1S t ne exchange interaction between the i — th 
Mn atom from the sublattice a and the Au atom, 
and x is t ne local spin susceptibility of the Au atom. 
Strong SOC on Au yields the MAE contributions due 
to M|^. For the uniaxial term it can be written 
as, E*jf a — —k^(M^ d ' z ) 2 , where the renormalized 
anisotropy constant fc|l is proportional to the square 
of the SOC parameter £ 2 d . Using Eq. ([2]) we can write, 
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Summing Eq. @ over the Mn sublattices a, we get 
zero magnetic moment which complies with the over- 
all time-reversal symmetry of the AFM band struc- 



K 4 ±,K Ml MAE constants (total and element specific in ture The total MAE due to Au, £o=t \E b / i 



meV per f.u.) of M^Au. Additional uniaxial MAE K 2 \i 
(in meV per f.u. per 1% strain) is induced by in-plane 
strain along the [100] crystal direction (see text). 

The magneto-crystalline anisotropy energy (MAE) 
is evaluated using the torque method which is imple- 
mented in the FP-LAPW basis [12J . For the tetrago- 
nal symmetry case, the phenomenological total energy 
dependence on the spin quantization direction reads, 
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non-zero, however, as seen from Eq. (|3j). It originates 
from a combination of strong 5c? SOC (fefl)) strong 
exchange splitting induced by the 3<i magnetic ele- 
ment ( J3d-5dM 3<i ), and the enhancement of the local 
spin susceptibility \- Eq. ([3]) implies that MAE is 
predominantly governed by two-site anisotropy and is 
proportional to the number of bonds between the Au 
and Mn atoms. 

Large magneto-crystalline anisotropics in bimetallic 
AFMs have already been reported in ab initio studies 
of the common AFM Mnlr. The key role of polarized 
orbitals of the noble metal has not, however, been 
identified in these works. To test that the physics 
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described in the previous paragraph is generic, we 
have examined the element specific MAE of Mnlr. 
In the calculations we considered the experimental 
lattice constants of the Ll Mnlr and the unit cell 
with two f.u.'s in a checkerboard collinear AFM struc- 
ture [II, [ll]- For the total MAE (K 2 ± + K AX + K A \\) 
we obtained -3.365 meV/f.u. (-6.73 meV per unit 
cell). This number is in good agreement with pre- 



vious results (-7.05 meV per unit cell in Ref. 1 19. -6.81 
meV per unit cell in Ref. 17). Our element-specific 
decomposition of the MAE shows that, similarly to 
the case of M^Au, the MAE is dominated by the 
contribution from the noble metal which again car- 
ries zero net moment in Mnlr. For the leading uni- 
axial MAE term K%± we obtained a contribution of 
-3.38 meV/f.u. from Ir and a much weaker contribu- 
tion of 0.024 meV/f.u from Mn. 

We now turn to the discussion of the magneto- 
transport effects in AFMs. The AMR effects are even 
in magnetic moment and, similar to the MAE, can 
be finite in the compensated AFM. In what follows 
we provide estimates of the AMR in the Coulomb 
blockade and tunneling regimes. To date TAMR and 
CBAMR have been observed only in fcrromagnets. 
We argue that both these effects are present and can 
be large in compensated AFMs. 

The CBAMR builds on a general principle that 
transport through an electronic device depends on po- 
sitions of electro-chemical potentials in the relevant 
electrodes. Transport and band structure with SOC 
are in this case directly linked via the anisotropy of 
the chemical potential with respect to the orientation 
of magnetic moments [HI EH- Single electron tran- 
sistors (SETs) are arguably the most sensitive de- 
vices to detect this effect. Large magneto- Coulomb 
oscillations of the conductance are induced by rotat- 
ing magnetic moments, whenever the changes in the 
chemical potential become comparable to the single- 
electron charging energy of the central island in the 
SET 10]. Our calculated difference between chemi- 
cal potentials for the staggered magnetization aligned 
along the [110] axis and the [001] axis is -2.5 meV 
for Mn2Au and 3.2 meV for Mnlr. This is compara- 
ble to chemical potential anisotropies in ferromagnetic 
semiconductor (Ga,Mn)As and in bimetallic ferromag- 
ncts 10]. The AFM CBAMR or a magnetoresistance 
of other AFM devices which are sensitive to changes of 
electro-chemical potentials of the order of a few mV's 
should therefore be readily detectable. 

For FM tunneling devices it has been demonstrated 
that a useful qualitative analysis of the differential 
TAMR can be obtained by considering its propor- 
tionality to the energy dependent anisotropy in the 
density of states (ADOS) in the magnetic electrode 
with respect to the crystallographic orientation of 
magnetic moments Q- Anisotropy in the group- 



velocity weighted tunneling density of states (AT- 
DOS) is considered instead of the ADOS when refer- 
ring to tunnel junctions with high quality interfaces, 
i.e., high degree of in-plane momentum conservation 
during tunneling [jj @] We emphasize that in either 
case these calculations are only approximate as they 
neglect the anisotropy of tunneling matrix elements 
between wavefunctions in the magnetic and nonmag- 
netic electrodes of the TAMR structure [l9| . 

Calculated energy-dependent ADOS and ATDOS 
for Mn 2 Au and staggered moments aligned with the 
easy [110] and hard [001] staggered magnetization 
axes are plotted in Fig. 2(a). The ADOS can be as 
high as 50% in the vicinity of the Fermi level and 
oscillates strongly with energy. This implies a sizable 
TAMR whose magnitude and sign are bias dependent, 
similar to the experimentally observed TAMR charac- 
teristics in bimetallic ferromagnets [H 0] . The calcu- 
lated ATDOS is also large and shows weaker depen- 
dence on energy around the Fermi level. In Mnlr the 
ADOS between the easy [100] and hard [001] staggered 
magnetization axes is weaker but still reaching 10%, 
as shown in Fig. 2(a). Note that for staggered mo- 
ment reorientations between the easy and hard axes, 
Mnlr is expected to display weaker TAMR of the two 
considered AFMs while the corresponding magnetic 
anisotropy constant K21. and the anisotropy in the 
chemical potential are stronger in Mnlr. This illus- 
trates that despite their common SOC origin the dif- 
ferent magnetic anisotropy phenomena can behave to 
some extent independently. 

The absence of a net magnetization makes AFMs 
rigid to magnetic fields and therefore one needs to re- 
sort to manipulation of the staggered moments via 
internal fields. One approach, previously reported in 
the literature, is based on exchange-spring effects [2(| 
occurring when the AFM is interfaced with a ferro- 
magnet. This technique can be used to trigger the re- 
orientation between the easy [110]-axis and hard [001]- 
axes which produces the TAMR effects predicted in 
Fig. 2(a). Instead of this exchange field method we 
now return to the magneto-crystalline anisotropy and 
demonstrate that the staggered moments can be con- 
trolled via this internal field without involving ferro- 
magnets. Specifically we consider the dependence of 
the MAE and of the corresponding in-plane easy-axis 
orientation in the AFM on the lattice strain. (This 
can be controlled externally by, e.g., a piezoelectric 
stressor f2lj|.) 

As shown in Tab. I, the in-plane easy-axis of un- 
strained M112AU lies along the [110] (or [110]) crys- 
tal direction. By applying a sufficiently strong strain 
along one of the cube edges, the easy-axis, and 
therefore also the staggered moments, will rotate to- 
wards the [100] or [010] direction, depending on the 
sign of the strain. To model this effect we elon- 
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FIG. 2: (a) DOS anisotropies for the hard ([001]) and 
easy ([110]) axes of M112A11, and the hard ([001]) and easy 
([100]) axes of Mnlr. (b) In-plane DOS anisotropies for 
staggered moment aligned along the easy axis ([010]) of 
the strained Mn2Au crystal and along the easy axis ([110]) 
of the unstrained Mn2Au. The energy is measured from 
the Fermi level. 



gated/contracted the unit cell along the [100]/[010]- 
axis, keeping the unit cell volume fixed. The corre- 
sponding crystal structure becomes orthorombic and 
the second order uniaxial anisotropy terms are given 
by sin 2 9[K 2 ± + K* cos 20]. 

Assuming strains up to 1% we found negligible 
changes in the out-of-plane uniaxial constant K 2 ± 
and in the fourth order constants K±± and -fQn- For 
the strain-induced in-plane uniaxial MAE we obtained 
K%n = 0.07 meV/f.u. per 1% strain and again the 
dominant contribution comes from the Au-atom (0.08 
meV/f.u. per 1%). The easy-axis shifts from the 
in-plane diagonal to the [010] direction when K 2 \\ > 
2-fQii . It means that only a fraction of a per cent strain 
is required to rotate the staggered moments between 
the [110] and [010] axes. These crystallographic di- 
rections are non-equivalent in the tetragonal Mn^Au 
and can yield a sizable TAMR, as shown by the cor- 
responding in-plane ADOS in Fig. 2(b). Note that 
the detailed TAMR characteristics can be modified by 
the strain as illustrated by the other ADOS curve in 



Fig. 2(b) which corresponds to the relative difference 
between the DOS in the unstrained case and moments 
along the [110] axis, and the DOS in 1% strained case 
and moments along the [010] axis. Nevertheless, the 
primary role of the strain is in the MAE where it trig- 
gers the reorientation of the magnetic easy-axis. 

In conclusion, we have proposed a set of relativis- 
tic magnetoresistance effects which can open a route 
towards AFM spintronics and which can be realized 
in AFMs alone without involving ferromagnetic el- 
ements. Our previous works in which TAMR and 
CBAMR were predicted for transition metal ferromag- 
nets 0, fToj ) using similar relativistic ab initio tech- 
niques have been subsequently confirmed by exper- 
iments [H, 0, [ll|. This together with the identified 
microscopic physics behind the magnetic anisotropy 
effects in compensated AFMs gives us confidence that 
predictions presented in this paper are realistic. 
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